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VO . ABSTRACT 

m ■ 

^^ , Recent results from the Planetary Nebula Spectrograph (PN.S) survey have re- 

■^ ■ vealed a rapidly falling velocity dispersion profile in the nearby elliptical galaxy NGC 

^^ ' 3379, casting doubts on whether this intermediate-luminosity galaxy has the kind 

OO , of dark matter halo expected in ACDM cosmology. Wc present a detailed dynamical 

^^ ■ study of this galaxy, combining ground based long-slit spectroscopy, integral-field data 

^ \ from the SAURON instrument, and PN.S data reaching to more than seven effective 

•'^ i radii. 

^^ ' We construct dynamical models with the flexible x^-made-to-measure particle 

H \ method implemented in the NMAGIC code. We fit spherical and axisymmetric mod- 

. - - i els to the photometric and combined kinematic data, in a sequence of gravitational 

potentials whose circular velocity curves at large radii vary between a near-Keplerian 
decline and the nearly flat shapes generated by massive halos. 

Assuming spherical symmetry we find that the data are consistent both with near- 
isotropic systems dominated by the stellar mass, and with models in massive halos 
with strongly radially anisotropic outer parts (/? > 0.8 at 7 Re). Formal likelihood 
limits would exclude (at la) the model with stars only, as well as halo models with 
Vcirci^Re) ^ 250kms~ . A sequence of more realistic axisymmetric models of different 
inclinations and a small number of triaxial tests confirm the spherical results. All 
valid models fitting all the data are dynamically stable over Gyrs, including the most 
anisotropic ones. 

Overall the kinematic data for NGC 3379 out to 7i?e are consistent with a range 
of mass distributions in this galaxy. NGC 3379 may well have a dark matter halo 
as predicted by recent merger models within ACDM cosmology, provided its outer 
envelope is strongly radially anisotropic. 

Key words: galaxies: elliptical and lenticular - galaxies: kinematics and dynamics - 
galaxies: individual: NGC 3379 - methods: iV-particle simulation - methods: numerical 
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1 INTRODUCTION 

There is strong evidence that most galaxies are surrounded 
by massive dark matter (DM) halos. This is most evident in 
spiral galaxies, where the rotation curves of extended cold 
gas disks remain flat out to large radii. In elliptical galaxies 
the evidence for dark halos has built up more slowly, and 
their halo properties are not so well known, because of a 
lack of ubiquitous tracer similar to the HI rotation curves in 
spirals. Only in a few cases is it possible to measure masses 
from extended HI ri ng velocities {e.g. iFranx et al.l 1 1994 : 
lOosterloo et al.ll2002l ). 

However, at least for giant elliptical galaxies 
stellar-dyn amical studies fro m int egrated light spec- 
tra {e.g. iKronawitter et all [2OO0I: iGerhard et al.l I2OOII : 



ICappellari et al.l I2OO6I : iThomas et al.l [20071), analyses of 



the X -r ay emitting hot ga s atmos ph eres (e.g. Awaki et al.l 
1994 [Matsushita et al.l 1 19981 : iLoewenstein fc Whitel 



I999I : iHumphrev et all I2OO6I : iFukazawa et all |2006| ) . and 



gravitational lens ing dat a 
Treu fc KoopmanJ |2004|: 



Koopmans et al.l 



(e.g. I Wilson et al.l 
Rusin fc Kochanekl 



20061 : iGavazzi et al.l 1200/1 ') 



200 ll: 



20051 : 



are now 



giving a fairly consistent picture. The general result from 
these studies is that these ellipticals are surrounded by 
dark matter halos, the inferred mass profiles (luminous 
plus dark) are nearly isothermal, i.e., the circular velocity 
curves approximately flat, and the dark matter contributes 
~ 10 - 50% of the mass within R^. The central DM 
densities in ellipticals are higher than in spir als, presumably 
reflecting their ear lier formation epochs (| Gerhard et al.l 
l200ll : lThoma3l2006l ). 

In light of this, the fin ding of lRomanowskv et al.l (|2003l ) 
and iDouglas et al.l (|2007l ) , that several intermediate lumi- 
nosity ellipticals (NGC 3379, NGC 4494, NGC 821) appar- 
ently have only difi'use dark matter halos if any, is quite 
surprising. Could the dark matter properties of these el- 
lipticals be different f rom those of giant ellipticals (e.g., 
iNapolitano et al.l |2005| ). perhaps related to the fact that 
these lower-luminosity galaxies are less often found in 
groups or clusters? The result of lRomanowskv et al.l (|2003l ) 
is based on the outer velocity dispersion profiles of the 
three galaxies, determined from individual planetary neb- 
ulae (PNe ) velocities measure d with the special PN.S in- 
strument (jDouglas et al.l [2003) • Two of the three galaxies 
are nearly round on the sky, and therefore the dynamical 
analysis was carried out with spherical models. A fourth 
galaxy with a fairly rapidly d eclining outer velocity dis- 
persion profile is NGC 4697 (JMendez et al.l bOOlh. How- 
ever, using axisymmetric particle models Ide Lorenzi et al.l 
(|2008l ) have recently shown that only models with massive 
dark halos are consistent with all the kinematic data for 
this galaxy, and that the best models have circular velocity 
Uc(5i?e) — 250kms~^ at 5 effective radii. Unfortunately, the 
diffuse gas envelopes of these intermediate luminosity ellipti- 
cals have very low densities, so an independent confirmation 
with X-ray data is difficult. 

The results of Romanowskv et al.l (|2003l ) were criticized 
bv lDekel et al.l (|2005l ). These authors pointed out that the 
well known mass-anisotropy degeneracy in the study of ve- 
locity dispersion profiles does not allow one to unambigu- 
ously determine the mass profile, that the triaxial nature of 
elliptical galaxies can cause low line-of-sight velocity disper- 



sions at some viewing angles, or that the PNe could trace 
young stars generated during the merger for mation instead 
of the bulk of the ol d stars as usually assumed . iDouglas et al.l 
(j2007i) argued that lRomanowskv et al.l (|2003l ) properly took 
into account orbital anisotropics in the data fitting process, 
that the effect of triaxiality is very unlikely to be present in 
all three galaxies, that the PN number density and veloc- 
ity dispersion profiles match the corresponding integrated 
light profiles reasonably well, and that this as well as the 
universality of the bright end of the PN luminosity function 
rules out that PNe only trace a young stellar population. 
IDouglas et al.l (|2007r ) concluded that their results continue 
to conflict with the presence of dark matter halos as pre- 
dicted in cosmological merger simulations. 

The issue is important enough to merit a further care- 
ful analysis. In this paper we construct dynamical mod- 
els of NGC 3379 with the very flexible NMAGIC particle 
code, making use of a variety of kinematic data, including 
SAURON integral-fleld data, slit kinematics, and the PN 
dispersion profile. 

The NMAGIC method is flexible not only with regard 
to anisotropy, but also in allowing axisymmetric or triax- 
ial shapes with radially varying axis ratios. This is impor- 
tant since the in t rinsic shap e of NGC 3379 is still in doubt . 
ICapaccioh et all (|l99J ) and lStatler fc Smecker-Hand (|l999l ) 
have suggested that NGC 33 79 is a triaxial SO galaxy seen 
almost face-on. IStatlerl (|2001j) considered triaxial dynamical 
models and constrained the shape of this galaxy to be ax- 
isymmetric and oblate in the inn er parts and triaxial in the 
outer parts. IShapiro et al.l (|200y) argue that the most likely 
model is one of a moderately inclined oblate system. 

The outline of the paper is as follows. In Section [2] we 
describe briefly how the various observational data for NGC 
3379 are used in the modelling. In Section [3] we give a few 
detaUs of the x^M2M NMAGIC method, and show how it 
performs on a mock galaxy data set similar to that for NGC 
3379. In Section|4]we then construct various dynamical mod- 
els for the real galaxy data, both spherical and flattened, in 
a sequence of potentials with increasing circular velocity at 
large radii. As summarized in the final Section [5] of the pa- 
per, our main conclusion is that the combined kinematic 
data for NGC3379 is consistent with a range of dar k matter 
halos including those found by iDekel et al.l (|2005|) in their 
cosmology-based merger simulations. 



2 OBSERVATIONAL DATA 

We begin by describing the observational data used in this 
study, which are all taken from the literature. We also de- 
scribe here the procedure employed for obtaining the three- 
dimensional luminosity density from the surface brightness 
data. In th e following we ad opt a distance 9.8 Mpc to 
NGC 3379 jjensenet al.ll2003l ). effective radius R^ = 47" 
(2.23 kpc), and an absolute B magnitude Mb = —19.8 
( Douglas et~aLll2007h . 



2.1 Photometric Data 



The photometric data used in the prese nt work consists 
of the wide-field B-band photometry of ICapaccioli et al.l 
(|l990l ). combined with the HST V-band observations of 
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Figure 1. C o mbin ed photometry of NG C 3379 from 
ICapaccioli et alj lll990l ) (open black circles) and ICebhardt et alj 
||2000|) (full red circles). The two panels show the surface bright- 
ness (SB) profile and the ellipticity e as a function of major axis 
distance. Beyond 500" the SB points arc from a Sersic model 
fitted to the interior data, and outside 81", the ellipticity has 
been set to e = 0.14. In the ellipticity panel the error bar with 
size 0.02 illustrates the typical errors in the outer ellipticity mea- 
surements. The isophotal shape parameters 04 and ag are not 
measured; they are set to zero. The lines show three-dimensional 
luminosity models determined from these data and reprojected 
onto the sky, for assumed inclinations of i = 90° (edge-on, full 
lines), i = 50° (dashed lines), and i = 40° (dash-dotted lines). 



iGebhardt et al.l ((20001) to increase the spatial resolution 
within the inner 10". The photometry has been matched up 
by assuming a constant color offset B — V = 1.03. The last 
eight surface brightness (SB) points from ICapaccioli et al.l 
()l99d ). outside 7? ~ 500", show fluctuations of an amplitude 
which we jud g ed un physical; these points we have replaced 
with a ISersid (|l968l ) proflle fitted to the galaxy further in. 
The same Sersic fit is used to extrapolate the SB proflle 
outside the last measured point at _R = 676". Similarly, we 
have replaced the measured ellipticities for R > 81", where 
the observational uncertainties become large, by e = 0.14. 
Figure [l] presents the combined photometric data, showing 
surface brightness and ellipticity e. The isophotal shape pa- 
rameters 04 and as are not available for these data and are 
thus set to zero. For the spherical models, we have used the 
SB proflle rescaled to a mean radius Rm = ^/ab — a\/l — e. 
For the axisymmetric models, we ha ve used a constan t PA o f 
70°; the isophotal PA measured by ICapaccioli et al.l (|199G| ) 
are within ±3° of this value. 



2.2 Deprojection 

In our implementation of NMAGIC a particle model can 
be fitted to the surface brightness and/or the deprojected 
luminosity density, cf. Section[3l Below we use both options, 



so first need to construct models for the three-dimensional 
luminosity density, j. 

In the spherical case the surface brightness can be de- 
projected uniquely. For an axisymmetric system the depro- 
jection is unique only for edge-on galaxies; for systems in- 
clined at an angle i with respect to the line-of-sight, the SB 
map contains information about the luminosity density only 
outside a "cone of ignorance" in Fourier sp ace, of openin g 
angle 90° — i, when i = 90° denotes edge-on (|Rvbickilll987l ). 
Thus, the deprojection of a moderately inclined galaxy re- 
sults in undetermined konus densitie s ([Gerhard fc BinnevI 
1 19961 : iRomanowskv fc Kochane^l 19971 ) . 

We deproject the surface brightness of NGC 3379, 
without correctin g for PSF effects, using the program of 
iMagorrianI (| 19991 ). The program finds a smooth axisymmet- 
ric density distribution consistent with the SB distribution 
for the specified inclination angle, by imposing that the solu- 
tion maximizes a penalized likelihood. This ensures that the 
shape of the 3D luminosity density is smooth and biases the 
model towards a specified diskyness. Because of the disklike 
nature of the undetermined konus densities, requesting the 
luminosity density to have zero diskyness effectively chooses 
between the different density distributions that fit the SB 
data for i 7^ 90° . We have used the program to compute lu- 
minosity densities for NGC 3379 for the inclinations i = 90° , 
i = 50° and i — 40°. Figure [l] compares the observed pho- 
tometry with the three deprojections reprojected onto the 
sky. Varying the inclination from 90° to 40° changes the 
intrinsic shape of the galaxy from El to E3. 



2.3 Kinematic Data 

2.3.1 Stellar-absorption line data 

We have taken long-slit absorpti on line kinematics from 
the li terature. We use data from IStatler fc Smecker-Hanj 
(|l999l ) at four different position angles, extending out to 
radii of ~ 80". We comp lement these kinematics with 
the spectroscopic data from iKronawitter et al.l (l200d) . The 
major axis slits fr om Statler fc Smecker-Hand (|l999l ) and 
IKronawitter et al.l (|2000l ) are misaligned by 10° in PA; 
however, the data along both major axis slits follow each 
other closely. The meas urements along the shifted slit of 
IKronawitter et al.l (|2000l ) reach 100" from the center. From 
both kinematic data sets we have the line-of-sight velocity, 
velocity dispersion, and higher order Gauss-Hermite mo- 
ment s hs and /14 (JGerhardl Il993l : Ivan der Marel fc Franxl 
I1993I ). Figure [5] shows the schematic arrangement of the 
kinematic slits used in the dynamical modeling. 

In addition to the long-slit kinematics we also use the 
integral-field spectroscopy obtained with the SAURON in- 
strument. Th e se kin ematic data were kindly provided by 
IShapiro et al.l (|200g) and consist of line-of-sight velocity, ve- 
locity dispersion and higher order Gauss-Hermite moments 
up to he. The SAURON field-of-view (FoV), shown by the 
(blue) rectangle in Figure (2] extends from —19.6" to 24.4" 
along its short boundary and from —34.8" to 35.6" along 
the long boundary. In this FOV, the positions of the 55x88 
"lenslets" with which spectra were taken define a fine grid 
of 4840 grid cells, which serve as the basis grid to define the 
1602 voronoi cells on which the final kinematic measure- 
ments are given. This results in a total of 9612 kinematic 
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Figure 2. Schematic view of tiie positions with kinematic data 
as used to construct the dy nam ical models. The slits f rom 
IStatler fc Smecker-Hanj l ll999l ) and lKronawitter et all (|2000| ) are 
coded in red and black, respectively. Boxes along the slits show 
the region of the galaxy for which respective kinematic data points 
were derived; these boxes are used to determine the luminosity- 
weighted Gauss-Hermite moments. The blue rectangle indicates 
the SAURON field-of-vicw. The ellipse shown is oriented along 
PA=70°, the average major axis of the photometry, and has a 
semi-major axis of length i?e and axis ratio q = 0.9. 



SAURON observables, as well as 1602 bin-luminosity ob- 
servables. The SAURON data are reproduced and compared 
to dynamical models in Section U) Each of the six panels 
shows the 1602 voronoi bins, giving (from left to right) v, 
(J, /i3, /i4, hs a nd he- A comparison of t he SAURON data 
with the data of lKronawitter et al.l (|200d ) along their major 
axis is given in Figure |3] Overall, the two data sets agree 
well with each other. The same is true for the compariso n 
of the SAU RON data withlstatler fc Smecker-Hand (|l999D . 
as shown bv lShapiro et al.l (200a). 

Both the SAURON data and the slit data are 
slightly asymmetric with respect to the center of 
the galaxy. If we denote the original SAURON 



I{x,y\vioa, iTios, hs, hi, /i5, he), 



and 



dataset with 

with /*(a;,i/|iiios,o-ios,/i3,/i4,/i5,/i6) = I{-x,-y,\ - 
viosjO'ioB,— ha, h4,—h5, he) the dataset obtained from I by 
point-symmetrical reflection with respect to the origin, we 
can construct a symmetrized dateiset I = 0.5(7 + /*). This 
symmetrized dataset I has a x^ per data point with respect 
to I of x^ /N = 1.01 when the original errors are used. Any 
point-symmetric model fit (spherical, axisymmetric, triax- 
ial) to the original data I will therefore have a systematic 
error floor of this magnitude. In the models below, we will 
actually fit the symmetrized SAURON data to avoid any 
systematic effects, but keep the o riginal errors on bot h sides 
of the galaxy separately (see also IShapiro et al.ll2006l ) . 

In a similar fashion, we have constructed symmetrized 
slit data sets. To do this we average the two points at nearly 
similar radius on both sides of the slit with respect to the 
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Figure 3. Comparison of the line-of-sight velociy distribution 
data along the galaxy's major axis {PA = 70°). The black circles 
correspond to t he SAURON data and th e open square symbols 
in red show the iKronawitter et al.l KOOOl ) data. The upper panel 
compares the original datasets, the lower panel is for the sym- 
metrized data. In each panel from top to bottom are shown: v, 
cr, /i3, h4, h^, and he, for the latter two there are only SAURON 
data. 



center. Taking into account the sign reversals of v and /13, 
we take for the symmetrized data point the weighted mean 
of the points on both sides, with weights proportional to 
the inverse square of the measurement errors, and assign a 
new weighted error for the averaged point. If a+ and a^ 
are the errors on both sides, the weights are w^ — l/cr+, 



ui_ = 1/(T_, and the new error a is given by the maximum 
of 2/(T^ = 1/(7^ + l/a^ and half of the deviation between the 
original data points on both sides. Again, the symmetrized 
data have a x^ jN — 1.0 systematic deviation from the orig- 
inal data, and therefore we will fit the symmetrized data be- 
low to avoid the model being pulled around by points with 
small error bars but large systematic deviations. The second 
panel of Figure [3] com pares the symmetrized SA URON data 
with the svmmetrized IKronawitter et al.l (|2000l ) data along 
the same slit as before. Again, the two data sets agree well 
with each other. 



2.3. 2 PNe data 

Planetary nebulae (PNe) are dying low- to intermediate 
mass stars that emit most of their light in a few narrow lines 
of which the [OIIIJA5007 is the most prominent one. Because 
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there are hardly other emission sources in elhptical galaxies, 
they can be detected fairly easily, and once identified, their 
line-of-sight velocity can be estimated from the Doppler shift 
of the emission line. The PN population in elliptical galaxies 
is expected to arise from the underlying galactic population 
of old stars and hence the PNe can be used as kinematic 
tracers for the stellar distribution. Their number relative to 
the luminosity of the galaxy is param etrized by the a pa- 
ram eter, which is a function of colour (JHui et al.lll995l ). 

iDouglas et al.l (|2007l ) processed observations of NGC 
3379 conducted with the Planetary Nebula Spectrograph 
(PN.S) instrument and detected 214 spatially and spec- 
trally unresolved PN candidates of which 191 are assigned 
to NGC 3379 . Usin g the "friendless" algorithm applied by 
iMerrett et al.l (|2003l ) they identified a small number of veloc- 
ity outliers, probably unresolved background galaxy contam- 
inants, which would be uniformly spread in velocity. The al- 
gorithm determined that 2 emission objects were more than 
n = 5 standard deviations away from the centroid of the 
velocity distribution of their A'^ = 15 nearest neighbours, 
an d 3 objects more tha n 3 standard deviations a (see Fig. 8 
of lDouglas et al.l 120071 ). The 3a line itself has considerable 
uncertainty at large radii, due to the small number of PNe 
found there. Thus the exclusion of the outermost outlier is 
somewhat uncertain. Because this object does have some in- 
fiuence on the outermost velocity dispersion point, we will 
compare the models to the data obtained both with and 
without this PN. 

The radial distribution of the PNe in NGC 3379 was 
found to be consistent with the stellar density profile, and 
their kinematics consistent with absorption-line data in the 
region where the data sets overlap. Because the kinematics 
of the PNe in NGC 3379 are dominated by random motions 
with little azimuthal variation, the velocity dispersion can 
be computed in radial annuli without losing significant dy- 
namical structure. We will thus use the radial run of the az- 
imuthally averaged PN velocity dispersion in the dynamical 
modelling, but also compare the models to the individual 
velocities in a relative likelihood sense (cf. the tables and 
figures in Section |3}. 



Gebhardt et all l2003l: iThomas et all |2004|: IValluri et all 



20041: ICappellari et al.1 l2006l : iChaname et all 12007 
van den Bosch et al.l |2008| ) the particle approach is rela- 



tively new and there are as yet only a few gal actic dynamics 
studie s in which it has been employed. iBissantz et al.l 
(1200J) mad e a first practic al app lication of the M2M 
method of ISver fc Tremaind (|l996l ) and constructed a 
dynamical model of the Milky Way's barred bulge and disk 
by constraining the projected density map. First attempts 
to extend the M2M method to account for kinematic 
observables in addition to de nsity constraints were made b y 
Ide Lorenzi et al.l (120061 ) and Ijourdeuil fc Emseiieinl (|2007h . 
However, a proper treatment of observat ional errors was 
not ye t included in their implementations. Ide Lorenzi et al.l 
(12007|) incorporated this in their x^M!2M algorithm and 
demonstrated the potential of the NMAGIC code by 
constructing particle models for spherical, axisymmetric, 
triaxial and rotating target stellar systems. Some exten- 
sions of the method and the first detailed modelling of slit 
kinematic and PN data for an elliptical galaxy (NGC 4697) 
are described in Ide Lorenzi et al.l (|2008l 'l. 

The NMAGIC method is fiexible not only with regard 
to the orbit structure, but also in allowing axisymmetric 
or triaxial shapes with varying axis ratios. Contrary to 
Schwarzschild's method, the "best" stellar density and lumi- 
nous potential need not be specified beforehand, but can be 
found from the evolution of the model. This makes it ideal 
for the present study because different intrinsic shapes have 
been suggested for NGC 3379 (see the Introduction) and the 
issue of whether the kinematics require or allow dark matter 
may well be connected not only with the orbital anisotropies 
but also with the detailed shape of the stellar density distri- 
bution of the galaxy. Given that NGC 3379 is nearly round 
on the sky, we have constrained most models in this paper 
to be axisymmetric, with density distribution fixed from the 
deprojection; however, for some models (in Sections 14.2.11 
and 14. 3|) the stellar system evolves towards a final, "best" 
density distribution, allowing for radial variations in axis 
ratio. This approach has proved sufficient for answering our 
main science question. 



3.1 Luminous and dark mass distributions 



3 NMAGIC MODELLING 

To investigate the amount of dark matter consistent with 
the kinematic data for NGC 3379, we construct a range of 
dynamical models for the stellar component of this galaxy. 
We use the fiexible x^-made-to-measure (x^M2M) particle 
me thod as described and implem ented in the NMAGIC code 



20081). Y^M2M is a deve lopment 



Sver fc Tremaind (|l996l ) that is 



bv lde Lorenzi et al.1 (|2007l . 
of the M2M algorithm of 
suitable for modelling observational data. The M2M meth- 
ods work by gradually adjusting individual particle weights 
as the model evolves, until the N-particle system reproduces 
a set of target constraints. In x^M2M the standard x^ statis- 
tics is used in the function to be maximized upon conver- 
gence of the weights. This allows for a proper treatment of 
observational errors, and the quality of the final model can 
be assessed directly from the target data. 

Compared to the familiar Schwa r zschild method (e.g. 
Schwarzschildl Il979| : iRix et al.l Il997l: Ivan der Marel et al.1 



19981 : ICretton et al.|[l999l : iRomanowskv fc Kochanekll20o1 



As in Ide Lorenzi et al.l (|2008f ). we assume that the luminous 
mass of NGC 3379 follows the light and characterize it by 
a constant mass-to-light ratio T, so that the stellar mass 
density is given by p^, — Tj. The total gravitational poten- 
tial is generated by the combined luminous mass and dark 
matter distributions, (j) = (j)* ~\~ 4'l, where (j>* is generated 
by p* = Tj. Only the luminosity density j is represented 
by the A^'-particle system. Its potential is compute d using a 
spheric al harmonic decomposit ion as described in ISellwoodI 
1:200^: lde~Lorenzi et all |20o3). The stellar potential is al- 
lowed to vary during the modeling process, but the DM halo 
is rigid. 

Here our aim is not to determine the detailed shape 
of the dark matter halo in NGC 3379, but rather to first 
see whether the PN velocities allow or require any dark 
matter at all in this galaxy. To answer this question we 
will investigate a one-dimensional sequence of potentials 
whose circular velocity curves vary at large radii between 
the near-Keplerian decline expected when the mass in stars 
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dominates, and the n early flat shapes gener ated by mas- 
sive dark halos. As in lde Lorenzi et alj ( 20081 ) we thus rep- 
resent the dark matter h alo by the logarithmic potential 
iBinnev fc Tremaine|[l987[ ) 



Mr) 



|ln(r^ + r^) 



(1) 



3.2 Model and target observables 

Target observables include surface or volume densities and 
line-of-sight kinematics. For modelling the luminosity distri- 
bution of NGC 3379, we generally use the deprojected lumi- 
nosity density of NGC 3379, expanded in spherical harmonic 
coefficients Aim on a 1-D radial mesh of radii r^. The corre- 
sponding model observables are computed from the particles 
based on a cloud- in-cell (CIC) scheme; see lde Lorenzi et alJ 
(|2007h . 

In some models, we do not constrain the three- 
dimensional luminosity density but only the stellar surface 
density, leaving the former free to evolve. In the remaining 
cases, we constrain the model by both the deprojected lumi- 
nosity density and the projected surface density. In a similar 
spirit as for the volume density, we use as target constraints 
for the observed SB distribution the coefficients of a Fourier 
expansion in the azimuthal angle, computed on a 1-D radial 
mesh of projected radii R^. For the corresponding model ob- 
servables, the particles are assigned to the radial grid using 
a CIC scheme, and the Fourier coefficients a™ and &,„ for 
the particle model on shell k are computed via 



Om,fc = L2_^lki cos{mifii)wi 



= L yj 7^/'^ sin(m<^i)u;i, m > 



(2) 



(3) 



where Wi are the particle weights, tpi their angular positions, 
and 7fe/ is a radial selection function. We use units for 
which the light Li of a stellar particle can be written as 
Li — Lwi with L the total luminosity of the galaxy. 

As kinematic constraints, we use the luminosity- 
weighted Gauss-Hermite coefficients from the SAURON or 
slit data, and luminosity-weighted v elocity moments fo r 
the PN data. For the SAURON data (|Shapiro et al.ll2006l ). 
the luminosity-weighted coefficients are determined from 
the truncated Gauss-Hermite representation of the line-of- 
sight velocity distribution (LOSVD) up to order /ig and 
the lumi nosity in the corresponding Voronoi bin. For the 
slit d ata (jStatler fc Smecker-Hanelll999l : iKronawitter et al.l 
l200d ) , they are constructed again from the measured Gauss- 
Hermite moments, up to order hi, and the luminosity in the 
slit se ction corresponding to the relevant LOSVD. The PN 
data (jDouglas et al.ll2007l ) are modelled either as 1-D radial 
dispersion profile or as a discrete set of velocities; in the for- 
mer case we use as suitable observables the second velocity 
moments vf^^, luminosity- weighted by the number of PNe 
per radial bin. 

The corresponding model observables j/j are construced 
from the particles via equations of the form 



Viii) = ^WiKj [zi(f)], 



where Wi are the particle weights and Zi are the phase- 
space coordinates of the particles, i — 1,---,N. Here the 
Kernel Kj corresponds to the observable yj. Detailed ex- 
pressions for the kinematic model observables are given in 
Ide Lorenzi et al.l (|2007l . |2008| ). 

In general, we replace the observables by the cor- 
responding temporally smoothed quantities to increase 
the effective number o f particles in the sy stem, cf. 
ISver fc Tremaind (119961 ): Ide Lorenzi et al.l ((2003). For the 
parameters chosen, the smoothing is typically over ~ 10'^ 
correction time steps. 



3.3 Constructing a particle model for the target 
data 

Generating an NMAGIC model for a set of observational 
constraints proceeds by evolving the force-of-change (FOC) 
equations for the particle weights. 



dwi{t) 
dt 



ewi{t) 



dS__M-^K, [z,(i)] 



^{Yi. 



A,(i) 



(5) 



(4) 



depending on the discrepancies between model (j/j) and tar- 
get observables (Yj), Aj(t) = {yj — Yj)/a{Yj). Here o-(Yj) 
in the denominator is the error in the target observable. 
Evolving the particle weights to convergence in this way is 
equivalent to maximizing the merit function 

F^fiS-^x' (6) 

with respect to the particle weights Wi, where for the profit 
function S we use the entropy, and the standard x^ mea- 
sures the goodness of the fit. The parameter fj, controls the 
contribution of the entropy function to F. The entropy term 
pushes the particle weights to remain close to their priors, 
so models with large /i will have smoother distribution func- 
tions than those with small p. The best choice for ^ depends 
on the observational data to be modeled, e.g. spatial cover- 
age and phase-space structure of the galaxy under consid- 
eration, but also on the initial conditions, and will be de- 
termined for the NGC 3379 dataset in the following Section 

EH 

Any NMAGIC model starts from a suitable initial 
model. For the model s presented in this paper, we have used 
as initial conditions a lHernguistl (|l990i ) model particle real- 
ization generated from a distribution function (DF), using 
the method described in lDebattista fc SellwoodI (|2000l ') . The 
particle realization consists of 7.5 x 10^ particles, has a scale 
length a = 1, maximum radius rmax = 60, and a total lu- 
minosity of unity. In model units, the gravitational constant 
is G = 1. In real units, the model lengthscale corresponds 
to 50". Thus when we match the model to NGC 3379, the 
effective radius of NGC 3379 becomes 0.94 model units, or 
2.23kpc at a distance of 9.8Mpc. 



3.4 Anisotropic mock galaxy model 

To prepare for the modeling of NGC 3379, we now con- 
struct a spherical mock galaxy model with known intrinsic 
properties to determine the optimal value of the entropy 
"smoothing" para meter /i in equa t ion ^ . Fo Uowing a simi- 
lar approach as in lCerhard et al.l ([1993) and lThomas et al.l 
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Figure 4. A comparison of the surface brightness profile of the 
mock galaxy model (dotted line) with that of NGC 3379 (full 
line), along the major axis of NGC 3379. 



(J2005f l we determine for which value of fi the fitted particle 
model best reproduces the intrinsic velocity moments of the 
input mock galaxy model. The "best" value of fi depends on 
the observational data to be modelled and their spatial cov- 
erage, on the phase-space structure of the galaxy, but also 
on the initial conditions from which the NMAGIC modelling 
starts. The same value can then be used for the modelling 
of NGC 3379, provided the mock galaxy is a reasonable ap- 
proximation to the real galaxy. 

For the lumin osity density of the mock galaxy we use 
a lHernguistl (|l990l ) model with total luminosity L — 1.24. x 
10^" Lq^b and scale radius a = 0.8 kpc, corresponding to 
Re « 30" for the distance of NGC 3379. As FigureHshows, 
the surface brightness profile of this model galaxy is a good 
approximation for NGC 3379. 

For the distribution funct ion of the mock galaxy we 
take an Osipkov-Merritt model (|Qsipkovll979l : lMerrittlll985l ') 
with anisotropy radius Va = 9a, giving an anisotropy pro- 
file similar to some of our later models fo r NGC 3379. 
The L OSVD kinematics is calculated following [Carollo_et_alJ 
()l995l ) and setting the mass-to-light ratio to Ts — 5. To the 
final LOSVD parameters we add Gaussian random variates 
with la dispersions equal to the respective error bars of the 
corresponding NGC 3379 measurements at that point. In 
this way we compute v, a, /13 and /14 points for the mock 
galaxy along all slits shown in Figure [2] Figure [5] compares 
the kinematics of NGC 3379 with the mock galaxy model 
along the major axis. 

In addition, we construct SAURON mock kinematics for 
each voronoi cell in the NGC 3379 data as follows. We first 
compute the velocity profiles as above at a few nearby radial 
positions. Using the spherical symmetry, we interpolate v, a 
and the higher order moments to the mid-cell positions of the 
fine grid described in Section[2]3] using a spline interpolation 
scheme. Then we compute the mock data for each voronoi 
bin by a luminosity weighted average over those cells of the 
fine grid which contribute to the voronoi cell under consid- 
eration. Finally, we add Gaussian random variates to the 




Figure 5. Comparison of the LOSVD kinematics of the mock 
galaxy model with those of NGC 3379 along its major axis, and 
with the particle model fit for /x = 2 X 10^. The red open s quare s 
show the NGC 3379 data from IStatler fc Smecker-Hand |l99i), 
the black circles show the pseudo data, and the solid line the self- 
consistent particle model obtained from fitting the pseudo data. 
The model data points are averages over the slit cells (see Fig.[2]l, 
and are connected by straight line segments. The panels from top 
to bottom are for v, cr, hg and /i4. 



kinematics with la dispersions corresponding to the respec- 
tive SAURON error bars in this voronoi bin. The SAURON 
pseudo data are shown in the top panels of Figure |6] 

The kinematic data set is completed with mock PN dis- 
persion data, using the projected velocity dispersions from 
ICarollo et al.l (|l995l ) and the errors from the observed PN 
dispersion points in NGC 3379. For test purposes, we also 
use two other sets of mock PN data with smaller errors (see 
Figure O. Again, Gaussian random variates corresponding 
to these errors are added to account for the scatter in the 
velocity dispersion points. Finally, we complete the mock ob- 
servational data set with the photometric constraints. In the 
entropy tests here, we restrict ourselves to spherical models, 
so in the expansion of the luminosity density the only non- 
zero term in the spherical harmonics series (c/. Section [3. 2|l 
is the radial fight in shells, Lk = V4^Aoo,k- However, to 
ensure sphericity, we also need to use the higher order coef- 
ficients A2o,k, ■ ■ ■, ^22, fc and Age,*! as constraints, set to zero. 
We define these photometric observables on a grid of radii 
rfe, quasi-logarithmically spaced in radius with inner and 
outer boundaries at rmin = 0.01" and rmax = 2500". We as- 
sume Poisson errors for the radial light a{Lk) = ^j L-kLIN 
where N is the total number of particles used in the par- 
ticle model and L is the total fight of the system. To esti- 
mate the errors in the higher order luminosity moments, we 
use Monte-Carlo experiments in which we compute the Aim 
many times from random rotations of a particle realization 
of the target density distribution. In these experiments the 
number of particles is 7.5 x 10^, which is the same number 
as in the x^M2M models. 

We now construct self-consistent particle models for the 
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Figure 6. Top panel: SAURON mock kinematic data for an anisotropic spherical galaxy model. Bottom panel: Self-consistent particle 
realization obtained from a model fit with /.t = 2 X 10'^. From left to right: u, a and the higher order moments h^-ha. 



anisotropic model galaxy target in a two step process, using 
the mock observations as constraints for NMAGIC. First, 
we start with the particle model described in Section 13.31 
and evolve it using NMAGIC to generate a self-consistent 
particle realization with the desired luminosity distribution 
(mock particle model), fitting only the photometric con- 
straints. Then, we use the mock particle model as initial 
conditions to fit both the kinematic and photometric target 
constraints for different values of n. 

The results are presented in Figure [S] The lower panel 
shows the goodness of the fit as a function of n, both in terms 
of the normalized x^ per data point and in terms of the merit 
function F from equation ([6}. The upper panel shows the 
rms relative difference A between the true internal velocity 
moments of the mock galaxy and those of the particle model 
realizations obtained for different values of fi. The intrinsic 
kinematics of the particle models are computed by binning 
the particles in spherical polar coordinates, using a quasi- 
logarithmic grid with 21 radial shells bounded by rmin = 
0.01" and rmax = 500.0", 12 bins in azimuthal angle (p, and 
21 bins equally spaced in sin 6. As can be seen from the top 
panel of Figure |8] there is no minimum in the rms A as a 
function of /i, but the particle models recover the internal 
moments of the input model well for /.i < 2 x 10'^. For larger 
H, the rms A increases rapidly because of oversmoothing 
in the model. The lower panel of the figure shows that x^ 
per data point is below unity for a large range of ^ but 
then increases for fi > 2 x 10^. In our modelling of NGC 



3379 below we have confirmed that this value of /i allows 
the models to converge towards strongly anisotropic orbit 
distributions. We have therefore used /x = 2 x lO'^ in Section 
[4] throughout. This is indicated by the solid symbol in Fig.|8l 

Figures [Sl|5] and [5] compare the anisotropic mock galaxy 
model and the particle model obtained from the target data 
with /i = 2 X 10'^. Figure \5\ shows the target kinematics 
along the galactic major axis and the corresponding par- 
ticle model kinematics. Figure [6] compares the SAURON 
mock data with the two-dimensional kinematics obtained 
from the particle model. Figure [7] shows the PN velocity 
dispersion data and compares to the model dispersion pro- 
files. All kinematic data are fit very well by the model. In 
fact, it is evident from Fig. [S] that the model is smoother 
than the mock data themselves, which is a consequence of 
entropy-smoothing and time-smoothing. 

Figure [^ shows how well the internal kinematics of the 
particle model for /i = 2 x 10'' compare with the intrinsic 
kinematics of the mock galaxy target. The velocity disper- 
sions (Tr, (T0 and ag, the streaming rotation v^, and also the 
anisotropy parameter f3g = 1 — al/a^ are reproduced well 
by the model within ~ IRe- At larger radii, the anisotropy 
of the Osipkov-Merritt distribution function cannot be en- 
tirely recovered even with idealized data (many dispersion 
points with small error bars), because there is no constraint 
from the data on the outer particle model (beyond the 
dashed line in Fig. [9] ). This is consistent with similar tests in 
iThomas et al.l (|200J). With fewer mock PN dispersion data 
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Figure 7. Different moek PN veloeity dispersion data for the 
anisotropic spherical target galaxy model, and corresponding 
model velocity dispersion profiles obtained with fi = 2 X 10"^. 
The three mock data sets differ by the quality with which they 
represent the underlying target velocity dispersion profile. The 
black dia monds have be en obtained using the original errors of 
the Dou glas et alj J2007|) PNe velocity dispersion data. The blue 
squares show the data obtained using reduced errors of 5kms^ , 
and the red circles illustrate an idealized data set with an in- 
creased radial resolution. The corresponding NMAGIC models 
are presented by the black dash-dotted line, the blue dash-dot- 
dot-dot line and the red dashed line, respectively. The latter is 
indistinguishable from the target velocity dispersion profile. 



points and larger errors, the particle model obtained from 
the target data and isotropic initial conditions is less tightly 
constrained; it becomes even less radially anisotropic despite 
fitting the actual PN data points well. 

Mass-to-light ratio So far all model fits have been made 
with the mass-to-light ratio fixed to the actual value used for 
the mock galaxy, T = 5. Now we investigate how accurately 
we can recover T with the dynamical models, given the spa- 
tial extent and quality of the observational data. To this end 
we fit particle models to the mock galaxy observations for 
different mass-to- light ratios in the range T £ [3, 10], keep- 
ing the entropy parameter fixed sd, fi = 2 x 10^. The results 
are presented in Figure [Till which shows how the quality of 
the model fit varies as a function of T, both in terms of 
X^ per data point and merit function F. As expected, the 



best model is obtained for T 
approximately 0.7. 



5; it has x per data point 



4 DYNAMICAL MODELS OF NGC 3379 

In this section we construct dynamical models for NGC 
3379 to learn about its stellar and dark matter distribution. 
We investigate mainly spherical and axisymmetric models, 
with and without dark matter halos, and fit the photome- 
try, SAURON integral field data, slit kinematics, and PNe 
velocity data. Our aim in this paper is not to constrain the 
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Figure 8. Top: Deviation rms A{/i) between the internal ve- 
locity moments of the final mock galaxy particle model and the 
input model. The upper and lower curves show the rms A(/i) 
obtained, respectively, with the original mock PN velocity dis- 
persion points, and with the errors of these dispersions reduced 
to Skms"'^. Bottom: The circles show x^ P^r data point of the 
model fit to the kinematic and photometric targets as a function 
of entropy parameter fi. The triangles display the merit function 
F, cf. equation ^. For both quantities, the curves obtained with 
the two sets of errors fall on top of each other. The full symbol 
indicates the optimal value chosen for fi. 
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Figure 9. Internal kinematics of the anisotropic target galaxy 
and its final particle model realizations. From top to bottom: ar, 
cr^f,, erg, ii0 and anisotropy parameter /3g = 1 — a^ /ai . The kine- 
matic quantities of the input mock galaxy are shown by circles, 
and are compared to those of three different particle models gen- 
erated for fi = 2 X 10'^: for idealized PN data (many points with 
small errors, dashed line), for mock data equivalent to the cor- 
responding NGC 3379 data (dot-dashed), and to NGC 3379-likc 
PN data with errors and scatter in the dispersion points reduced 
to 5kms~^; see Fig. [7] 
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Figure 10. Recovering the mass-to-light ratio of the mock galaxy. 
The quality of the model fit as a function of mass-to-light T is 
shown in terms of x^ P^r data point (circles) and merit function F 
(triangles) . All models are generated from the target pseudo data 
for fi = 2 X 10^. The input mass-to-light ratio T = 5 is recovered 
as the minimum in the curve x^/^(1f")i where the model fit has 
xVN ~ 0.7. 



detailed halo mass profile of the galaxy, but only to ascer- 
tain whether a dark matter hal o is allowed, or required , by 
the kinematic data. Thus, as in lde Lorenzi et al.l (J2008l ) we 
investigate a simple sequence of potentials which include the 
self-consistent part from the stellar component and a fixed 
halo potential as in equation (ITJ. The circular speed curves 
corresponding to these potentials vary at large radii from 
the near-Keplerian decline expected when the mass in stars 
dominates, to the nearly flat shapes generated by massive 
halos. They are shown in Figure [TT] and their halo potential 
parameters are given in Table [l] 

In the following subsections, we describe spherical mod- 
els (*J4T| and oblate models ( i^4.2|l . as well as a few models 
without imposed axisymmetry constraints ( i]4.3|) . and then 
discuss the significance of the fits to the data in a separate 
subsection ( i^4.4|l . To begin with we construct self-consistent 
particle models for NGC 3379 in which the distribution 
of stars is spherica l. This allows for an easy comparison 
with previous work (JRomanowskv et al.ll2003l : lDouglas et al.l 
I2OO7I ). 



4.1 Spherical models 

4. 1.1 Target data and modeling process 

First we must determine the photometric and kinematic 
observables. Analogous to Section 13.41 we use the spheri- 
cal harmonics expansion coefficients Ai^^ of the deprojected 
luminosity density as target data to constrain the particle 
models. Specifically, we use j4oo, A20, A22, • • •, ^66, but set 
all terms higher than Aoo to zero, adopting the same radial 
grid as in Section 13.41 Errors for the luminosity terms are 
estimated as in Section 13.41 As kinematic observables, we 
use the SAURON and slit kinematics, as well as the binned 
FN velocity dispersion profile; see Sections 12.31 and 13.21 The 
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Figure 11. Circular velocity curves for the potentials used in 
the dynamical modelling, including the self-consistent stars-only 
model A (dashed line), and models including different spherical 
dark matter halos in addition to the stellar component (solid 
lines, from bottom to top: models B, C, D and E). For this figure 
the distribution of stars is assumed to be spherical with mass- 
to-light ratio as given by the final NMAGIC model for the data 
in the respective spherical potential. The shaded area shows the 
ra nge of circular velo city curves in the merger models discussed 
bv lDekel et"al] |20oi); see Section [5] 



SAURON data and most slit data are symmetrized, only the 
slit parallel to their minor axis of lKronawitter et al.l (|2000l ) 
cannot be symmetrized and for this slit the original kine- 
matic data points are used. 

We match the particle models to these data in the fol- 
lowing three-step process, (i) We begin with the initial par- 
ticle distribution described in Section [33] and evolve it with 
NMAGIC to a self-consistent model that reproduces the tar- 
get Aim. (ii) Starting with this density model we then con- 
struct dynamical models, fitting the full set of photometric 
and kinematic target observables. If the potential includes a 
dark matter halo, we first relax the density model for 1000 
steps in the total gravitational potential (c/. Section l3.ip . 
assuming a mass-to-light ratio of 8, to make sure that the 
model is in approximate equilibrium. After this relaxation 
phase, we evolve the particle system for ~ 10^ NMAGIC cor- 
rection steps while applying the complete set of constraints. 
During the correction phase the mass-to-light ratio T is ad- 
justed i n parallel, us ing its own force-of-change equation as 
given in lde Lorenzi c t al. (2003)- After each correction step, 
the potential generated by the particles is updated but the 
dark matter potential (if present) is constant in time. In this 
process the entropy parameter has value /j, — 2x 10''; cf. Sec- 
tions [33]a-iid 2113] (iii) In the final step, we keep the global 
potential constant and evolve the system freely for another 
5000 steps, without changing the particle weights (phase- 
mixing). This completes the modeling process. Thereafter 
we generally evolve the model with all potential terms ac- 
tive for a further 10000 steps to test its stability. For refer- 
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0.641 
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2231.9 


8.03 
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3175.2 


6.73 



E* 



260 



0.241 



0.084 



0.215 



0.522 



0.504 



2649.4 



3.52 



Table 1. Table of parameters and fit results for models of NGC 3379 with spherical potentials. Models A-E correspond to the circular 
rotation curves in Fig. 1111 Model D+ is the same as D but for a higher value of the entropy. Model E* is the self-flattened oblate 
model in halo E of Section 14.2.11 For those models columns (l)-(3) give the model code and the parameters tq and vq used in equation 
JTJ for the respective dark halo potential. The next six columns list the x^ values per data point, for all obsorvablos [column (4)] and 
for the luminosity density and surface brightness constraints, the SAURON kinematic obscrvables, slit kinematic observables, and PN 
observables separately [columns (5)-(9)]. Column (10) gives the numerical value of the merit function in equation ^, and column (11) 
the final (B-band) mass-to-light ratio. The respective number of constraints are A' = 12997 for A-E and N = 12557 for E*, whith 
Nain, = 640, N,t, = 200, N,au = 11214, N,i = 1135, Npt^ = 8. 



15 



20 



25 



30 



35 




1000 



Figure 12. Comparison of the surface brightness profiles of the 
rcprojocted spherical models with the photometric data (points). 
The lines are for the spherical models A-E and the self-flattened 
model E*. 



ence, 10000 correction steps in the self-consistent potential 
correspond to « 110 circular rotation periods at R^, or 5.8 
Gyr. 



J^.l.Z Results 

In this way we obtain spherical dynamical models for NGC 
3379, reproducing the density expansion and all kinematic 
data including the PN velocity dispersion profile. Model A is 
the self-consistent model without dark matter halo, models 
B-E have halos of increasing circular velocities, as shown in 
Figure[TT]and Table[T] The quality of fit for these models can 
be judged from Table [l] which gives the numerical values of 
the merit function F and fists various values of x^ per data 
point, both those obtained globally for all the data, and 
those found for each of the four data sets separately (density 
expansion, SAURON, slit, and PNe). The 3a outlier point 



discussed in Section 12.3.21 is not included in the modeling 
and in the Xpn in Table[l] but its influence will be discussed 
below. 

Figures [T2ll 1 51 compare the different data with the mod- 
els. Fig. ll2l shows the surface brightness profiles, Fig. ll3l the 
integral field LOSVD parameter fields, Fig.[T3]the kinemat- 
ics along several slits, and Fig. [15] the PN velocity disper- 
sion profiles. The model SB profiles fit the observed profile 
very well, and agree with each other within the thickness 
of the lines in the plot. The SAURON data are fitted with 
xiau/Neau — 0.2 by all our spherical models. Notice that 
the particle noise in the models is significantly smaller than 
the noise in the symmetrized Sauron maps. Also the xii/^si 
for the combined slit data are less than unity; the plots for 
models (B,D) in Fig. 1141 show a few small systematic devia- 
tions but generally the fits are very good. In the central 30 
arcsec the slit data are dominated by the SAURON data. 
Notice that these spherical models are not constrained to 
be spherically symmetric also in their kinematic properties; 
hence they can also fit the observed (< 50kms^^) rotation 
of NGC 3379 with high accuracy. The small x'^-values are 
caused by the fact that the observational errors are slightly 
larger than the point-to-point fluctuations (see Fig. [3|, and 
to a greater extent, because with the entropy scheme we 
cannot smooth the models too much without erasing their 
anisotrop ic phase-space structure; cf. Sections 13. 41 14.1.^ and 
[dc Lorcn zi et al.l (|2008l ). 

The comparison of the models to the PN.S data is shown 
in Fig. [TE\ If we use the outermost dispersion point as given 
in lDouglas et al.l ([20071), models A-D with no or moderately 
massive halos provide a good match to the data, but the 
most massive halo model E fits less well, being high by ~ 
2cr with respect to the outermost dispersion point and by 
~ 1.3cr with respect to the second-outermost point. If we 
include the objec t classified as 3a "f riendless" outlier (see 
Section [2. 3. 2 1 and [Douglas et al.[[2007[ ) in the outermost bin, 
the Icr error range of the outermost PN dispersion point 
extends to significantly larger velocities; see the red open 
circle and error bar in Fig. 1151 Then model E also fits the 
PN dispersion profile, overestimating the outermost velocity 
dispersion point by less than la. 

The intrinsic kinematics of these spherical models 
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Figure 13. Symmetrized SAURON kinematic data for NGC 3379 (top row) compared with similar data extracted for the spherical 
models B and D and the self-flattened model E* (lower three rows). Notice that the particle noise in the model panels is significantly 
smaller than the noise in the corresponding data fields for all LOSVD parameters shown. In the panels for cr and /14 a slightly colder 
ring-like structure with larger /i4 hints at some deviations from spherical symmetry. 



is shown in Figure 1161 One recognizes the expected 
signature of the well-kn own mass-anisotropy degeneracy 
(jBinnev fc MamonI ligS^ ): In the more massive halos, the 
same falling line-of-sight dispersion profile requires larger 
radial anisotropy. Thus in the models with halo the ra- 
dial anisotropy rises outside l-2i?e- Particularly the more 
massive halo models D and E require strongly radially 
anisotropic orbit distributions (/3 ~ 0.9) to be consistent 
with the falling dispersion profile of NGC 3379. Radial 
anisotropy was suggested as one of the possible causes for 



the measured profile by iDekel et al.l (|2005l ). based on a 
comparison with their merger models. However, the typical 
anisotropics in their models are more moderate (/3 ~ 0.5). 

Despite their strong radial anisotropy, the massive halo 
models D and E show no sign of an instability when evolved 
freely after the model fitting and phase-mixing. Rather, they 
evolve very slowly, reaching after 5.8 Gyr of evolution a con- 
figuration with slightly triaxial shape (e < 0.1) in which the 
initial slow rotation has mostly gone away. A similar evo- 
lution is seen for the near-isotropic model A without dark 
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matter halo, indicating that this evolution may be connected 
to these equlibria being spherically symmetric only in their 
mass distribution but, due to the rotation, not in their kine- 
matics. In any case, the PN dispersion profiles do not change 
during the evolution, i.e., the constraints on the dark matter 
halo remain as before. 

In conclusion, the results of this section show that both 
near-isotropic spherical models with low density dark matter 
halos, and radially anisotropic spherical models with massive 
halos provide excellent fits to the available kinematic data 
for NGC 3379, including the PN dispersion profile to ~ 77?^. 
A more quantitative discussion is deferred to Section [4. 41 



4-1.3 Entropy smoothing 

The entropy term in the force-of-change equation ([5]) 
smoothes the particle models by trying to maintain the val- 
ues of the particle weights near their priors, here chosen 
as 1/N . Because all models start from an isotropic system 
with equal weight particles, the entropy smoothing thus bi- 
ases the final models towards isotropy and slow rotation. To 
allow the models to develop strong radial anisotropy in their 
outer parts requires a relatively low value of the entropy pa- 
rameter (see Section I3.4|l , which is below that appropriate 
for an isotropic system. Otherwise the constraints from the 
small number of PN dispersion points with their relatively 
large Poisson error bars would be overwhelmed by the en- 
tropy smoothing. We demonstrate this in Fig. [15] and Table 
[T]with a model D+ constructed with /j, = 2x 10*; this model 
is indeed degraded in its ability to fit the PNe data, relative 
to model D which is for the standard fj, — 2 x 10'^ in the 
same halo. 

Contrary to second derivative regularisation, say, en- 
tropy smoothing does not distinguish between local and 
global uniformity of the particle weights; it likes to have all 
particle weights similar to their priors. Thus if n is chosen 
such as to allow large difi^erences in weight between radial 
and circular orbits, it also allows similar differences between 
particles on neighbouring orbits if this is preferred by the 
data. With fi — 2x10 the models can therefore fit the data 
with x^/N < 1 as seen in Table [l] The efi'ect is strongest 
for the spherical models because these have a larger number 
of independent orbits than less symmetric systems. How- 
ever, Fig. [16] shows that the intrinsic velocity moments are 
smooth functions of radius, and below we will see that also 
the LOSVDs are smooth functions. Thus the good fits of 
the various models to the PN data are not achieved by large 
local variations of the orbital weights for orbits near the PN 
data points. 
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4.2 Oblate models including dark matter halos 

There is some evidence t hat N GC 3379 may be non- 
spherical. ICapaccioli et al.l (|l99ll ) argued that the bulge of 
NGC 3379 is remarkably similar to the one of NGC 3115, 
a well-known SO galaxy. Further, also the SAURON kine- 
matic data, shown in the upper panel of Figure 1131 show 
signatures of non-sphericity, particularly, a faint cold ring 
visible in the velocity dispersion and hi panels with pro- 



Figure 14. Comparison of models B (dashed lines), D (full lines) 
and E* (dash-dotted lines) with the symmetrized slit data along 
the major and minor axes from IStatler &: Smecker-Hanj (119991 ) 
(top and m iddle panel) and the unsy mmetrized minor-axis paral- 
lel slit from lKronawitter et al.l l(2000|) (bottom panel). The model 
data points are averages over the same slit cells as the target data 
(see Fig. I2}, and are connected by straight line segments. 
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Figure 15. Comparison of tlie PNe velocity dispersion profiles 
of the spherical mo dels with the PN.S data. The PNe velocity 
dispersion points of lDouglas et al.l 12003) ^'I'e shown as black di- 
amonds; when the object classified by them as 3(t "friendless" 
outlier is included, the la error range of the outermost PN dis- 
persion point extends to significantly larger velocities (red open 
circle). The dashed line shows the self-consistent particle model 
A. The solid lines represent the dynamical models including a 
DM halo, i.e., from bottom to top models B, C, D, and E. The 
heavy dotted line is for the higher-entropy model D+, and the 
dash-dotted line is for the self-flattened model E* . 




r/R, 



Figure 16. Intrinsic kinematics of the final spherical models A 
(dashed lines) and B,D,E (full lines), and the self-fiattened model 
E* (dash-dotted lines). Panels from top to bottom show the ra- 
dial, azimuthal, and vertical velocity dispersion profiles, the mean 
azimuthal streaming velocity, and the meridional anisotropy pro- 
file, all computed in an equatorial plane through the model (per- 
pendicular to the rotation axis for the spherical models, and per- 
pendicular to the line-of-sight for model E*). The models in the 
more massive dark matter halos are more radially anisotropic, as 
expected. The shaded area in the lower panel corresponds to the 
range of anisotropy profiles f ound in the el l iptica l galaxy remnants 
in the merger simulations of lDekel et al.l 1120051 ). 



jected radius R « 15'lj. Thus to understand how much dark 
mass around NGC 3379 is allowed by the kinematic data for 
this galaxy may require more general models than spheri- 
cal ones. In this section we will present oblate axisymmetric 
models in the family of spherical halo potentials considered 
already in the last section. 



4-2.1 Face-on oblate model in a spherical potential 

As a first step we attempt to construct a model for NGC 
3379 in a massive dark halo, in which the distribution of 
stars is flattened along the line-of-sight. This model is re- 
quired to have a small line-of-sight velocity dispersion at 
large radii, thus wil l be flattened in accordan ce with the 
virial theorem {e.g. iBinnev fc Tremainelll987^ . We do not 
know beforehand what the required shape of this model must 
be, so we will use the NMAGIC method to find it for us. 
Throughout this experiment the gravitational potential is 
constrained to remain spherically symmetric, being the sum 
of the spherical part of the luminous matter potential and 
the spherical halo potential. For illustration we embed this 
model in halo E, and will hence hereafter denote it as model 
E*. 

To construct this model we replace the Aim constraints 
(c/. Section 14.11) . which before imposed a spherical shape 



^ As can be seen from Fig. 1131 the feature can also be reproduced 
in spherical models. 




500 
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Figure 17. Surface brightness contours of the kinematically de- 
projected oblate model E*, when viewed perpendicular to the 
line-of-sight, i.e., edge-on. The model's outer parts have been pref- 
erentially fiattened to match the falling PN velocity dispersion 
profile. 
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on the particle distribution, by the Fourier moments of 
the surface brightness distribution given in Figure [l] They 
are computed from the photometry as in Section 13.21 on 
a grid in projected radius quasi-logarithmically spaced be- 
twen i?min = 0.01" and i?max = 1500". The higher-order 
moments are set to zero, enforcing axisymmetry. We then 
start from spherical initial conditions and use NMAGIC to 
flatten the particle model through fitting the kinematic ob- 
servables, particularly the PN velocity dispersion profile. As 
kinematic constraints, we use the SAURON, slit, and PNe 
velocity dispersion data. The entropy parameter is kept at 
the same value as for the spherical models, /x = 2 x 10'^ . Dur- 
ing this "kinematic deprojection", the spherically averaged 
potential generated by the particles is updated after regu- 
lar time intervals, but the non-spherical terms are ignored. 
The DM potential is given by equation ^ and remains con- 
stant in time. After the correction phase, the model is again 
allowed to freely evolve for some time. 

Figures [lH [H [13] and [15] show how the final "self- 
flattened" particle model E* compares to the various data. 
The model fits the data as well as the best-fitting spherical 
models. As anticipated, the model makes the PN dispersion 
profile compatible with a massive dark halo potential by 
flattening the outer distribution of stars and decreasing the 
model a along the line-of-sight. Fig. [16] shows that the line- 
of-sight velocity dispersion measured at ~ 2 — IRe in the 
equatorial plane [ae) is half the 0- dispersion in this plane; 
the radial dispersion still dominates, however. The model's 
flattening is illustrated in Figure 1171 which shows the SB 
distribution in an edge-on projection perpendicular to the 
line-of-sight. The axis ratio is q ~ 0.7, but to match the 
decreasing line-of-sight velocity dispersion profile, the flat- 
tening increases at large radii. 

While this model illustrates the power of the NMAGIC 
method, and provides an excellent fit to the photometric 
and kinematic data in a massive dark DM halo, it is not a 
realistic model for NGC 3379. For it is only in a spherical 
potential as assumed for model E* that a face-on distribu- 
tion of stars can show rotation. More realistic axisymmetric 
models must therefore be inclined to allow for the rotation 
seen in the SAURON and slit data. 



4.2.2 Self- consistent oblate models 

Therefore we now consider oblate models for NGC 3379 with 
inclinations i = 90°, i — 50° and i = 40°, in which the ax- 
isymmetric gravitational potential of the stellar component 
is computed self-consistently from the particles. We investi- 
gate models without DM as well as models including various 
DM halos as detailed in Table [2] The gravitational potential 
of the DM halo is still assumed to be spherical and is kept 
fixed. All models are evolved to fit the photometry, the slit 
and SAURON absorption line kinematic data, and the PN.S 
velocity dispersion profile. 

The procedure employed for constructing the models is 
similar to that in Section [4.11 Again, we need to first spec- 
ify the observables. We expand the deprojected luminosity 
density of NGC 3379 for each inclination in a spherical har- 
monics series and determine the expansion coefficients on 
the same quasi-logarithmic grid in radius as before. As ob- 
servables we use the luminosity on radial shells L^ and the 
higher order moments A20, A22, ■ ■ ■, Aee, but set the m 7^ 



terms of the expansion to zero to force the models to re- 
main axisymmetric. Errors for the Aim coefficients are esti- 
mated as in Section [3.41 We thus obtain three different sets 
of luminosity density observables Aim with corresponding 
errors, one for each of the three inclinations. In addition to 
the Aim, we also use the surface brightness itself as a con- 
straint, through the Fourier moment observables on the grid 
of projected radii Rk as in the previous Section 14.2.11 Er- 
rors for these Fourier moments are computed similarly as 
the Aim errors. The kinematic constraints are identical to 
those used for the spherical models: they are the luminosity 
weighted, symmetrized Gauss-Hermite moments from the 
slit data and SAURON data (see Section [2]3]), and the PNe 
kinematics represented by the binned line-of-sight velocity 
dispersion points. 

For the combined set of observables we construct par- 
ticle models in a similar three-step process as for the spher- 
ical models, (i) We start with the spherical particle model 
described in Section 13.21 and use NMAGIC to generate an 
equilibrium model with the desired luminosity distribution, 
as given by the deprojection of the photometry for the given 
inclination, (ii) We then use the resulting particle model as 
a starting point to generate the flnal set of models by fitting 
the photometric and kinematic constraints in the different 
DM halos. We use the same entropy parameter fj, = 2 x 10'^ 
as for the spherical models, (iii) Finally, we first keep the 
potential constant and let the system evolve freely without 
changing the particle weights, and thereafter test the stabil- 
ity of the model. 

The quality of the fit for the different halo models and 
inclinations is again characterized by the value of the merit 
function F of equation ([6]) and the values of the different 
X^ per data point, both globally and for the individual data 
sets. These are given in Table[2]and will be discussed further 
in Section [4.41 In addition to the models shown in Table [2] 
we have also constructed a similar suite of models for the un- 
symmetrized SAURON and slit data. These models were of 
similar quality as the models for the symmetrized data, i.e., 
when subtracting the systematic error floors determined in 
Section[2XI](xL^./iVsa^.(sys) = 1.0 and X%/Nsi{sys) = 1.0) 
from the x^ values of the models for the unsymmetrized 
data, the model x^ values became very similar to those re- 
ported in Table [1 

Figures I18I20I compare some of the final axisymmetric 
particle models to the SAURON, slit and PNe data. Both 
edge-on and inclined models again are very good matches 
to the SAURON and slit data, with or without dark matter 
halo. The PN velocity dispersion profile is fitted well by 
the models with the lower mass halo models B,C; halo D 
slightly overestimates the oute r PN velocity dispersion point 
given by [Douglas et al.l (|2007l ) but is within la of the outer 
point when the "friendless" outlier is included. Model E90 is 
incons istent with the outer dispersion point of iDouglas et al.l 
(|2007l ). cf. Table [5] but is only marginally inconsistent with 
the data when the outlier is included. Based on this together 
with the likelihood results reported below, halo D is the most 
massive halo consistent with the PN data. Figure \7T\ shows 
that for this model the dark halo contributes about 60% of 
the total mass within the radius of the last PN data point 
at ~ 7Re ~ 15kpc. 

Figure [22] shows the intrinsic velocity dispersions, 
streaming velocity, and anisotropy for some of the models. 
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Table 2. Tabic of parameters and x^~^^ results for oblate models of NGC 3379. Columns (l)-(3) give the model code and the parameters 
ro, vo used in equation JTJ for the respective dark halo potential; all halo potentials are spherical [q^ = 1.0). The fourth column gives the 
inclination i and the next six columns list the x^ values per data point, for all observables [column (5)], and for the density constraints, 
surface brightness constraints, SAURON kinematic observables, slit kinematic observables, and PN observables separately [columns (6)- 
(10)]. Column (11) gives the numerical value of the merit function in equation ^, and column (12) the final (B-band) inass-to-light 
ratio. The respective number of constraints are A^ = 13237, Naim = 680, N^t = 200, Nsau = 11214, A^^,; = 1135, Npn = 8. 



Because of the small projected ellipticity of NGC 3379, and 
the assumed spherical dark matter halo, the edge-on mod- 
els are very similar to the spherical models in the respec- 
tive halo potentials and the higher circular velocity halos 
require large radial anisotropy to match the PN data. The 
inclined flattened models have similarly small a^ — (re in 
the model equatorial plane, but somewhat larger a^, as ex- 
pected. Also in the axisymmetric models it is the radially 
increasing, strong radial anisotropy which causes the rapidly 
decreasing PN velocity dispersion profile in the massive dark 
halo potentials. 

Finally, we comment briefly on the stability of these 
models. All models in halos A-C show no signs of any change 
after 5.8 Gyr of evolution following the phase-mixing after 
the NMAGIC fit. The D models are almost unchanged, de- 
spite the strong radial anisotropy, developing after 5.8 Gyr 
a percent-level triaxiality just outside the error bars of the 
A22 constraints. The models in halo E show a similar slow 
evolution during which they in addition develop significant 
positive hi across the entire image. As in the spherical mod- 
els, the PN dispersion profiles remain unchanged during this 
evolution. 



4.3 Triaxial models 

We have constructed a small number of models for which 
the stellar density was not constrained to remain axisym- 
metric, in order to see whether the larger freedom in the 
orbit structure of non-axisymmetric potentials would allow 
the models to fit the PN kinematics also in the most massive 
halo E. However, these models do not have isophote twists: 
we have kept the constant value PA= 70deg for the posi- 



tion angle in the photometry, neglecting the observed small 
variations APA = ±3deg. These models are generated as 
follows, using the full power of NMAGIC: Only the surface 
brightness and kinematic data are used as constraints, in a 
similar way as for model E*, leaving all density Aim terms 
and corresponding potential terms free to change during the 
evolution. This allows the model to freely change its orienta- 
tion. As initial conditions we have used a spherical model, a 
model fiattened along the line-of-sight, or the inclined model 
D40. 

Because we know that valid models in halo D can be 
found, we first evolve a model in halo D, starting from ini- 
tial conditions D*, a model that had previously been ob- 
tained exactly analogously to model E* (see Section [4. 2. If) . 
Because of the line-of-sight streaming velocities, this system 
rotates out of the sky plane while NMAGIC simultaneously 
keeps adjusting the orbit structure to match both the sur- 
face brightness and the projected kinematics. This model 
converges to an almost axisymmetric model with inclina- 
tion i ~ 46 deg, and is then completely stable over 5.8 Gyr 
of evolution. We have computed iteratively the mean in- 
termediate and minor axis lengths inside ellipsoidal radius 
' + {y/bf + {z/cfY''^, following [Dubinski & Carlbcr3 



(|l99ll ). We obtain axis ratios b = 0.988 and c = 0.74 for 
q = 50" ~ 2.375kpc, and b = 0.990 and c = 0.73 for 
q = 200" ~ 9.5kpc. Note that the error in these axis ratios 
is about 0.002, due to the large number of particles used 
in the diagonalization of the tensor. This weakly triaxial 
model matches all the kinematic data, SAURON, slit, and 
PNe, very similar to models D40 and D50, and is listed in 
Table [2] as model DR. Its projected kinematics are shown in 
Figs. I18II20I and the enclosed dark matter fraction is about 
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Figure 18. Comparison of axisymmctric and weakly triaxial models with SAURON kinematic data for NGC 3379 (top panel). Following 
panels are for models A90, D90, DR. Model A90 has all the mass in the stars, while the later two models include a massive halo; see 
Table [2] 



60% of the total mass within the radius of the last PN data 
point at ~ 7i?e ~ 15kpc (Fig. [HI- 

Also shown in the Table are the results for model ER, 
which was obtained analogously starting from model E*. 
This model does not fit the PN data. None of our other 
attempts to obtain a valid model E has been successful, in- 
cluding one insp ired by some old work on merger remnants 
IjGerhardl Il983al lbf) , following which we tried to construct 
an oblate-triajcial model whose inner oblate parts are seen 
edge-on by the observer, while its triaxial outer regions are 
observed along the short axis. 



We believe the main reason for the failure in halo E 
is the observed rotation of NGC 3379, of which either the 
sense (along the projected major axis) or the amplitude do 
not allow the low-inclination configurations required by the 
low values of velocity dispersion at large radii. Consider a 
triaxial model viewed approximately along the short axis, 
which could easily accomodate the outer falling dispersion 
profile by a corresponding change of shape with radius as 
in model E*. Such a model is not consistent with the ob- 
served rotation because in this case the rotation visible to 
the observer would be around the model's long axis, i.e., it 



18 F. De Lorenzi et al. 



■-50 

250 

^200 

hso 

'100 

50 

0.2 

0,1 

" 

-0.1 

-0.2 

0,2 



-0.2 
-0.4 



^ 50 

; 

■-50 

250 

^200 

hso 

'100 

50 

0.2 

0,1 

" 

-0,1 

-0,2 

0,2 

:, 
-0,2 
-0.4 





-50 50 

R (arcsec) 



-^-^i-jJ= 



h— i-^- 



:-J-*^- 




-50 50 

R (arcsec) 

Figure 19. Comparison of the axisymmetric models A90 
(dashed Unes), D90 (full lines) and the weakly triaxial model 
PR (dash-dotted lines) with the symmetrized slit data from 
IStatler fc Smecker-Hang l|l999l ) along the major (top) and minor 
axes (bottom panel) . The model data points are averages over the 
same slit cells as the target data (see Fig. [J} ; and are connected 
by straight line segments. 



would be observed along the minor axis on the sky, whereas 
the actual observed rotation is along the major axis on the 
sky. On the other hand, a triaxial model with minor axis in 
the sky plane must correspond to a very round model unless 
it is near-prolate and viewed end-on. Then the box orbits 
are also viewed near end-on, making it difficult to arrange 
a falling outer dispersion profile. Thus valid triaxial models 
are likely to be radially anisotropic and to have inclinations 
near i ~ 45 deg as model DR, or be more edge-on and quite 
round, so the dominant effect on the inferred potential is 
the radial anisotropy. If so, the failure of model ER then 
suggests that extending our suite of models to include tri- 
axial halos (beyond the scope of this paper), is unlikely to 
increase the range of allowed circular velocities. 
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Figure 20. Comparison of the radial velocity dispersion profile 
from the PN.S data with the oblate and weakly triaxial particle 
models. The dashed line shows the stellar-mass only model A90. 
The other broken lines show models B90, C90, D90, the solid line 
shows model DR, and the upper dotted line shows model E90. 
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Figure 21. Enclosed DM fraction as function of radius for the 
final particle models B90, C90, D90, DR. 



4.4 Likelihoods and quality of the fits to the data 

We now turn to discussing the question which models are 
acceptable fits to the data and which models can be ruled 
out. To do this, it is customary to determine Ax'^ values 
relative to the best-fitting models, and determine the con- 
fidence boundaries according to the number of parameters 
to be determined. In our case, we essentially determine only 
one parameter, the halo circular velocity at ~ 7 Re, or vq, 
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Figure 22. Intrinsic kinematics of the final models A90, B90 
(dashed), C90, D90 (full) and D50 (dash-dotted lines). Panels 
from top to bottom show the radial, azimuthal, and 9 velocity 
dispersion profiles, the mean azimuthal streaming velocity, and 
the meridional anisotropy profile /3g = 1 — u'^/u'^. The models 
in the most massive halos are strongly radially anisotropic, as 
expected. 



SO the relevant Ax^ = 1 (the mass-to-light ratio of the 
models is optimized together with the weights). However, 
all our models match the Sauron and slit kinematic data 
to within la per data point, i.e., formally better than the 
underlyin g "true" model (cf. the discussion in Sections 13.41 
14.1.31 and Ide Lorenzi et al.l (|2008l )). Clearly, we cannot ap- 
ply a Ax^ = 1 for small variations within Icr relative to, 
say, the Sauron data points. Even if the best model fitted 
with exactly Xsau — 10*, this would make little sense: for 
Ngau = 10*, Ax'^ = 1 corresponds to an average change per 
data point of ~ lO"*^. Only if the Ax^ = 1 arises because 
of significant mismatch of a few crucial data points would 
this seem reasonable. The crucial data points for the issue 
addressed in this paper, the dark matter halo in NGC 3379, 
are the PN velocities or the binned PN dispersions. Thus we 
focus our discussion on the merit of the models relative to 
these data. 

Figure [53] shows the Xpjv S'UcI Ax^ values for both the 
spherical and the axisymmetric models from Tables [l] and (2] 
For the PN dispersion points we have 7 degrees of freedom 
(8 data points minus 1 fitted parameter), so expect y^ — 
8.18 (68.3% probability) for a typical good model. Thus we 
consider any model that fits the PN velocity dispersions to 
better than x^ ~ 8.18 as valid as the underlying "true" 
model and compute Ax^ relative to Xpjv — 8.18. The curves 
in Figure [23] are plotted for the two cases with and without 
the "friendless" outlier of lDouglas et al.l (|2007|) contributing 
to the outermost dispersion point. The models with halos A- 
D are allowed in both cases, while the models for halo E are 
consistent with the data only when the outlier is included. 

So far we have compared the models only to the PN 
velocity dispersion profile, rather than to the LOSVDs or 
unbinned velocity data. Figure [24] shows the LOSVD his- 
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Figure 23. x'^ S'^kI Ax^ values of the various spherical and ax- 
isymmetric models for NGC 3379 with respect to the PN velocity 
dispersion data. The full black lines (red dashed lines) connect 
the x^~v 3-lu6s obtained witho ut (including) the Scr "friendless" 
outlier of lDouglas et al.l 1120071 ). Ax'^ is computed relative to the 
expected value of x^ = 8.18 for 7 degrees of freedom. 
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Figure 24. Comparison of the PNe LOSVDs in the circular an- 
nuli corresponding to the outermost three bins in the velocity 
dispersion profile, with the LOSVDs of the nearly isotropic, low- 
density halo model B (dashed lines), and the radially anisotropic 
massive halo model D (solid lines), in the same circular annuli. 
The ordinate is in units of PN number, and the model LOSVDs 
have been normalized to the same integral over the velocity range 
shown. In the middle panel for the second annulus, the mean ve- 
locity of the PNe is non-zero at the ~ 3(t level; both models are 
inconsistent with this velocity distribution. The velocity distribu- 
tions in the other two panels are fully consistent with both mod- 
els. In the lower panel for the outermost shell, the PN histogram 
is shown with (dash-dotted) and without (solid l ine histogram) 
the 3c r-"' friendless"' outlier; see Section 12.3.21 and [Douglas et al.l 
||2007|) . 
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without outlier 


with outher 


Halo 


ln£ 


2Aln£ 


ln£ 


2Aln£ 


A 


-605.14 


2.50 


-611.18 


4.58 


B 


-604.21 


0.64 


-609.67 


1.56 


C 


-603.89 


0.0 


-608.89 


0.0 


D 


-604.74 


1.7 


-609.23 


0.68 


E 


-607.16 


6.54 


-611.38 


4.98 


A 


-608.50 


4.23 


-613.53 


5.01 


B 


-607.01 


1.25 


-611.82 


1.60 


C 


-606.38 


0.0 


-611.02 


0.0 


D 


-606.68 


0.60 


-611.14 


0.23 


E 


-608.81 


4.85 


-613.04 


4.02 



Table 3. LikeUhood values for the PN data in the spherical mod- 
els. Column (1): model code. Columns (2,3): log likelihood ln£ 
and difference 2Aln£ relative to the best model C, for the PN 
sample not includi ng the 3cr "friendless" outlier in the outermost 
shell, according to iDouelas et al.l 1 120071) . Columns (4,5): same, 
but for the PN sample including this outlier. The top half of the 
table refers to posterior likelihoods of the models fitted to the 
PN velocity dispersion profile, the lower half gives likelihoods for 
similar m odels in which the PNe were fitted with the likclhood 
method of lde Lorenzi et aiX l|2008l '). 



tograms for the PNe in the outermost three circular annuli 
used for computing the PN velocity dispersion profile, su- 
perposed on the LOSVDs of models B and D in the same 
radial shells. In the plot for the outermost bin, the PN his- 
togram and model LOSVD are shovy n with and without th e 
3(T "friendless" outlier according to IPouelas et al.l l)2007h . 
Both the near-isotropic low-density halo model B and the ra- 
dially anisotropic massive halo model D are consistent with 
the PN velocity distributions in the first and third annuli, 
and both appear inconsistent with the apparent non-zero 
mean motion of the observed PNe in the second annulus. 

Table |3] shows the posterior likelihoods of the spheri- 
cal models for the observed PN velocity data set, evaluated 
from the model LOSVDs in the eight radial shells used in 
the fits. Also listed are the likelihoods resulting from direct 
likelihood fits of the sph erical models to the PN data, using 
the method described in lde Lorenzi et al.l ([20081). Figure [25] 
shows a plot of these likelihoods as a function of the mod- 
els' circular velocity at 7i?e, the radius of the outermost 
PN dispersion point. Despite the small number of poten- 
tials investigated and the issue of whether the Scr "friend- 
less" outlier should be included, the overall shape of the 
likelihood function £ is not too far from the theoretically 
expected Gaussian. Thus we can determine a confidence in- 
terval from the condition A log C > 0.5, resulting in approx- 
imately 165kms~^ < Vcirc{7Re) ^ 250kms~^ at la. This 
would exclude both model A without dark matter and the 
most massive halo model E. However, we do not believe this 
is a very strong result, given the infiuence of a single outlier 
on the likelihood values in Table [3] and the asymmetries in 
some of the LOSVDs (see Fig. [21}. Note also that all mod- 
els are consistent with the data at the 2a level, for which 
VcirciTRe) <280kms-^ 
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Figure 25. Relative likelihoods from the data in Table [3| as a 
function of the model circular velocity at 7 Re. Open symbols 
show posterior likelihoods of models fitted to the binned disper- 
sion profile, full symbols show likelihood values based on direct 
likelihood fits to the PN velocities. Squares show likelihoods for 
the PN sample without t he Sa "friendle s s" ou tlier in the out- 
ermost shell, according to JDouglas et al.l 112007^ . circles for the 
sample including this outlier. 



5 SUMMARY AND CONCLUSIONS 

In this paper, we have carried out a dynamical study of the 
elliptical galaxy NGC 3379. This intermediate luminosity 
El galaxy has a rapidly declining velo city dispersion pro- 
file, w hich has been taken a s evid ence bv lRomanowskv et al.l 
(|2003l ') and [Douglas et all (|2007l ) that this galaxy may lack 
the kind of dark matter halo that the current ACDM cos- 
mology requires. 

To explore this issue further, we have combined pho- 
tometry, long slit spectroscopic data, SAURON absorption 
line kinematics and PN velocity dispersion data, to fit dy- 
namical models in a sequence of potentials whose circular 
velocity curves at large radii vary between a near-Keplerian 
decline and the nearly flat shapes generated by massive ha- 
los. The combined kinematic data set runs from the center 
of NGC 3379 to about 7 effective radu. 

For constructing the dynamical models we have used the 
flexible y^-made-to- measure particle code NMAGIC devel- 
oped by (|de Lorenzi et al.ll2007l . [2008h . The NMAGIC mod- 
els described in this paper consist of 7.5 x 10^ particles, and 
for the first time are constructed for such a comprehensive 
data set, including integral field kinematic data. 

We find that a variety of dynamical models both with 
and without dark matter produce viable fits to all the data. 
For assumed spherical symmetry we find that the data are 
consistent both with near-isotropic systems which are dom- 
inated by the stellar mass out to the last kinematic data 
points, and with models in massive halos whose outer parts 
are strongly radially anisotropic (/3(77?e) ^ 0.8). In these 
latter models, the stellar mass distribution dominates in the 
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center, and the dark matter fraction is ~ 60% of the total 
at 7 Re. 

In the spheri cal potentials we have also used the like- 
lihood scheme of Ide Lorenzi et alj (|2008l ) to fit the models 
directly to the PN velocities. From the likelihood values ob- 
tained in these fits as well as the posterior likelihoods of 
the models fit to the dispersion profiles, we estimate confi- 
dence limits on the halo circular velocity at IRg , resulting in 
approximately 165kms~^ < Vcirc{TRe) ^ 250kms~^ at la. 
This would exclude both the model without dark matter and 
the most massive halo model E in our sequence which has 
Vcirc{7Re) — 275kms~^. 

To illustrate the power of NMAGIC we have used it 
to find the shape of a model flattened along the line-of- 
sight in a spherical potential including this most massive 
halo E, which flts all the kinematic data with high accuracy. 
However, all attempts to find more realistic models with this 
massive halo have failed, suggesting that we may have found 
the upper limit of the range of consistent mass distributions. 

We have then constructed self-consistent axisymmetric 
models of inclinations i = 90° , i = 50° , and i — 40° in the 
same sequence of halos potentials. These models essentially 
confirm the spherical results. The edge-on models are very 
similar to the spherical models, becoming highly anisotropic 
in the more massive halos. The inclined models in addi- 
tion become more flattened at large radii, which helps in 
decreasing the outer velocity dispersion proflle. Finally, we 
have constructed a weakly triaxial model, in the most mas- 
sive allowed halo D, by free evolution from an axisymmetric 
model flattened along the line-of-sight. This model ended at 
i ~ 45 deg inclination, is almost axisymmetric, and matches 
all kinematic data very well, similarly to the inclined ax- 
isymmetric models in this halo potential. All these models 
are stable over Gyrs. 

Our main conclusions are as follows; 

(i) The kinematic data for NGC 3379 out to 7Re are 
consistent with a variety of potentials and do not give strong 
constraints on the mass distribution in this galaxy. The main 
reason for this is the well-known degeneracy between mass 
and radial anisotropy which is substantial when the velocity 
dispersion profile falls with radius. Formal confidence lim- 
its on the halo circular velocity at 7 Re are 165kms~^ < 
Vcirc{7Re) ^ 250kms^^ at la, which would weakly exclude 
models without dark matter. 

(ii) NGC 3379 may well have the kind of dark mat- 
ter halo consistent with the current ACDM paradigm. The 
cir cular velocity curv es of the merger models constructed 
bv iDekel et al.l (|2005h in the ACDM cosmology framework, 
Vcirc{r)/vcirc{Re) — {R/ Re)''''''^'^^ , fall right into the range 
of circular velocities of our best-fitting models in Fig. 1111 

(iii) Such models, however, are required by the data 
to have strongly radially anisotropic orbit distributions in 
their outer regions, /3 > 0.8 at 7Re, while mod el predictions 
at 7 Re are 3 < P <~0.6 JDekel et al.ll2005h and /? ~ 0.7 
IjAbadi et al.ll2006l ). Kinematic data at even larger radii than 
presently available would be required to discriminate be- 
tween these models and less anisotropic models with lower 
mass halos. 
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